Covariant Newtonian and Relativistic 
dynamics of (magneto)-elastic solid model 
for neutron star crust. 

B. Carter^ E. Chachoua^ N. Chamel^* 



t LuTh, Observatoire de Paris, 92195 Meudon, France, 
Copernicus Astronomical Center (CAMK), 00-716 Warsaw, Poland. 

1 July, 2005 



Abstract. This work develops the dynamics of perfectly elastic solid model for appli- 
cation to the outer crust of a magnetised neutron star. Particular attention is given to 
the Noether identities responsible for energy-momentum conservation, using a formu- 
lation that is fully covariant, not only (as is usual) in a fully relativistic treatment but 
also (sacrificing accuracy and elegance for economy of degrees of gravitational freedom) 
in the technically more complicated case of the Newtonian limit. The results are used 
to obtain explicit (relativistic and Newtonian) formulae for the propagation speeds of 
generalised (Alfven type) magneto-elastic perturbation modes. 



1. Introduction 

In astrophysical contexts of the kind exemplified by a neutron star crust, it is com- 
monly necessary to go back and forth between relativistic models having the advantage 
of greater elegance and in principle - particularly at a global level - of higher accu- 
racy, and Newtonian models that are more convenient from the point of view of other 
considerations such as economy in gravitational degrees of freedom, and availability of 
detailed underlying descriptions at a microscopic level. As a consequence of the fact 
that - unlike the Galilean invariance group - the Lorentz group is semi-simple, there 
are contexts (e.g. involving superfluidity [1, 2]) in which a fully relativistic treatment is 
actually easier to implement than a corresponding Newtonian treatment, even though 
the latter would be perfectly adequate as far as accuracy is concerned. On the other 
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hand there arc many contexts - particularly those involving electromagnetic effects or 
strong gravitational fields - in which a relativistic treatment would be indispensible if 
high accuracy were required, but for which a Newtonian treatment might nevertheless 
be easier to implement as a first approximation. 

To facilitate transition between these two complementary kinds of description, it is 
desirable to develop a unified treatment in which both relativistic and Newtonian models 
are described in terms of technical machinery and terminology that are compatible as 
much as possible, so as to be consistent in the limit when the relativistic spacetime 
metric goes over to the degenerate spacetime structure of Newtonian theory. In a 
coherent approach of this kind, the Newtonian limit will naturally be obtained in fully 
covariant formulation of the kind [3, 4, 5] whose mathematical machinery was first 
developed by Cartan [6]. In a preceding series of articles [7, 8, 9] on multiconstituent 
fiuid and superfiuid dynamics, it was shown how such a 4 dimensionally covariant 
formulation of Newtonian theory can provide physical insights (e.g concerning concepts 
such as helicity) that are not so easy to obtain by the traditional approach based on a 
3+1 space time decomposition. 

Continuing in the same spirit, the purpose of the present work is to contribute to 
the further development of the unified treatment of relativistic and Newtonian theory 
by treating the case of elastic solid models, of the kind appropriate for the description 
of the outer crust of a neutron star, including the magneto-elastic case (that arises in 
the limit of perfect electrical conductivity) for which the elastic structure is modified 
by a frozen-in magnetic field, of the kind whose effects are observed in pulsars. (The 
category of such models includes the special case of ordinary - fiuid not solid ~ perfect 
magnetohydrodynamics in the limit of negligible elastic rigidity.) 

Accurate description of such stars at a global level (not to mention a recently pro- 
posed cosmological application [10]) requires a general relativistic treatment, but use of 
a fiat space background will be sufficient for treatment of the local mechanical proper- 
ties to be considered here. In ordinary pulsars the magnetic field is sufficiently low that 
(except in the outer skin and the magnetosphere outside, where the matter density is 
comparatively low) such a flat background space treatment can be carried out (as de- 
scribed below) within a purely Newtonian framework. However a fully (at least special, 
if not general) relativistic treatment will be indispensible even locally (in a Minkowski 
background) when the magnetic field is sufficiently strong, as will be the case not just 
near the surface, but even in the deeper layers, for the special class of pulsars known as 
magnetars. The relativistic version of the magneto-elastic treatment developed here is 
particularly relevant for such strongly magnetised {B ^ lO^^G) neutron stars, in which 
fiares powered by magnetic stress are believed to be responsible for gamma ray bursts 
of the brief but intense kind observed in soft gamma repeaters, the most spectacular 
example so far - the most intense ever observed in our galaxy - having been the 27 
December event that occurred in SGR 1806 — 20 in 2004 [1 1] . For a complete treatment 
of such a flare a fully general relativistic description would presumably be necessary 
since a phenomenon of this kind is though to be attributable to a global modiflcation 
of the magnetic field of the neutron star [12, 13]. 

As an application, in both the Newtonian and fully relativistic cases such that the 
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underlying solid is in a simple isotropic state, the relevant (rigidity modified) propaga- 
tion speeds of (Alfven type) perturbation modes are explicitly evaluated. 

A subsequent article will be needed to treat the more elaborate kind of model 
needed for the innermost layers of a neutron star crust, in which an ionic sohd lattice 
is permeated by an independently moving current of superfluid neutrons. 

2. Milne structure of Newtonian spacetime. 

Before proceding, let us recapitulate the geometric essentials of Newtonian space 
time structure in a 4-dimensional background with respect to an arbitrary system of 
local coordinates x^, ^ = 0, 1, 2, 3, as described in greater detail in the first article [7] 
of the preceding series. 

The specification of a relativistic structure is fully determined by a non degenerate 
spacetime metric tensor having components 51^^, which in the special relativistic case 
are required to be constant in a prefered class of Minkowski type coordinate systems. 
This tensor will have a well defined contravariant inverse g^^, from which a correspond- 
ing connection (which will vanish with respect to Minkowski coordinates in special 
relativity) with components T^'^^ that are unambiguously obtainable using the usual 
Riemannian formula T^^^ = g'^'^{ga{fi,p) — ^/2gfip,a) using a comma to denote partial 
differentiation and round brackets to indicate index symmetrisation. 

Newtonian theory is traditionally formulated in terms of an Aristotelian frame, 
meaning a direct product of a 1-dimensional trajectory parametrised by the Newtonian 
time t and a fiat 3-dimensional Euclidean space whose positive definite metric gives 
rise to a corresponding 4-dimensional metric r;^^ that is of degenerate, rank-3, positive 
indefinite type, so that it does not determine a contravariant inverse tensor, and that, 
unlike its relativistic analogue, is not physically well defined because it depends on 
the choice of the Aristotelian "ether" frame, as characterised by the choice of a unit 
ether fiow vector, say, that will be a null eigenvector of the corresponding degenerate 
metric, i.e. that will satisfy r]^i,e^ = 0. There is however a complementary tensorial 
"Coriolis" structure that (unlike the Aristotelian structure) is physically well defined 
in the sense of being preserved by the allowable (time foliation preserving) ether gauge 
transformations. This invariant structure consists of the time gradient 1-form = 
and a contravariant metric tensor '-/^'^ that (like its gauge dependent covariant comple- 
ment r7^,^)has the property of being degenerate, of rank-3 positive indefinite type, with 
null eigendirection determined by the time covector, i.e. 

r''u = Q. (1) 

Although (like the non-degenerate metric in the relativistic case) this associated 
pair of tensors and 7'^^ is physically well defined, nevertheless the specification of 
this "Coriolis" structure (unhke that of the non-degenerate metric in the relativistic 
case) is not by itself sufficient to fully determine the geometric structure of spacetime 
in the Newtonian limit, and furthermore (like the metric in special relativity but not 
in general relativity) it is not freely variable over spacetime. 

To start with, the Coriolis structure specified by the fields 7'*'' and is required 
to be flat in the sense that there should exist coordinates with respect to which the 
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corresponding field components arc constant. It follows more particularly that there 
will be coordinates of a preferred type with respect to which these fields will have 
components of standard Aristotelian-Cartesian form, as given in terms of Kronecker 
notation by '^^^ = 5^5]^ + 5^5^ + ^zK-' ^^"^ ^ Such a coordinate system will 
determine a corresponding Aristotelian "ether" frame vector with components = 5^ 
and its associated covariant space metric with components r]^y = Sj^^Sl + S^S^ + Sj^Sf,, 
and = 5^. (In the preceding work [7] this quantity ?7^j, was written as 7^^, using the 
same Greek letter gamma as for its contravariant analogue, but in the present work the 
symbol 7 will be reserved for quantities that are gauge independent.) The Aristotelian- 
Cartesian coordinate system will also determine a corresponding symmetric connection, 
namely the one whose components 

r/p = V. (2) 

will vanish in that system, a requirement which evidently ensures that it will have 
vanishing curvature. The corresponding covariant differentiation operator V will satisfy 
the commutation relation 

VpV,-V,Vp = 0, (3) 

and will be such that the corresponding covariant derivatives of the tensor fields 
and t/j, will vanish: 

V,7'^'^ = V,t^ = 0. (4) 
As well as satisfying the algebraic conditions 

eH, = 1 e'^T},. = , rf^pY" -V;-S;- e% , (5) 

the ether frame dependent fields e'^ , 1]^^, and the associated space projection tensor 
r/^ introduced in (5), will also have corresponding covariant derivatives that vanish: 

Vp77^. = Vpr); = 0. (6) 



These fields can be used to specify an ether frame dependent Lorentz signature met- 
ric defined by = —tfjty + rjui,, with contravariant inverse, given by g'^'' = —e^e^ + 'j^'^, 
whose determinant provides a 4-dimensional volume measure that (modulo a sign ambi- 
guity depending on a choice of parity orientation) fixes a corresponding antisymmetric 
tensor with components e^^P^ — £^i^^p'^\ (using square brackets to denote index antisym- 
metrisation). This measure tensor is alternatively definable directly by the condition 
that its components with respect to Aristotelian-Cartesian coordinates (with the chosen 
orientation) should be given by or -1 whenever the indices are respectively even or 
odd permutations of the sequence {0, 1, 2, 3}. The corresponding antisymmetric covari- 
ant mcasTirc 4-form e^yp„ = ^[t^upcr] then specifiable (in the manner that is familiar in 
the context of relativistic theory) by the normalisation condition 

^"'^''^W = -4! ■ (7) 

It can be seen that (unlike and rj^j^y and the frame dependent Lorentz metric ) they 
share with and the property of being independent of the choice of ether gauge. 
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These tensors will evidently give rise to purely spacelike 3 index analogues, defined as 
the Hodge type duals of and e'^ respectively, namely the gauge independent space 
alternating tensor 

e""" = e'^'P'^t^ , (8) 

and the frame dependent 3-form 

^Hvp ^ S-Xnup (9) 

which is interpret able as representing an ether current. It is evident that all these 
tensors will share the flatness property characterised by the connection, in the sense of 
satisfying the conditions 

Va£^'^''^ = 0, Va£^.,<, = 0. (10) 

and 

V,e'^^p = 0, ^xe%,^0. (11) 



Despite the simplification provided by the flatness property that is expressed by 
(4), (6), (10) and (11), the Newtonian case is subject to the complication that neither 
the fiat coordinate system nor even the corresponding flat connection is unambigu- 
ously determined by the tensor fields 7'^'' and t^. The standard form expressible by 

^tiu _ _|_ ^ix^u _|_ _ ^jjj ^^^^ preserved by a large category of 

transformations that is known as the Coriolis group, which includes not only boosts but 
also time dependent rotations. However the physical structure of Newtonian spacetime 
is not preserved by time dependent rotations, but only by transformations of a more re- 
stricted category known as the Milne group, which is characterised by the condition that 
the modification of the ether frame vector should depend only on the Newtonian time t, 
as specified modulo a choice of origin by t^^ = 1^. This means that the transformation 
of the ether vector will be expressible in the form 

^e^ + , (12) 

for a boost vector field 6^ that is subject to the condition 

Ift^^Q, Y^^pb^ = (13) 

and for which the corresponding acceleration vector will be specified by 

at' = ePVplf, a% = 0. (14) 

The ensuing transformation of the covariant space metric will be given by 

Vt^u ^ Vf,^ = rip,y - 2ti^p,riy)pV + rip„¥h''tpty , (15) 

while those of the corresponding space projection and space measure tensors will be 

Tll^vl- b% , el^p ^ el^p = el^p + h^e^^^p . (16) 
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Unlike the ether vector e'', the covariant space metric rj^^, and the frame dependent 
Lorentz signature metric that was invoked above, the tensors constituting the Cori- 
ohs structure, namely the time gradient covector t^, the contravariant space metric j^^, 
and also the associated space-time measure given by e^i^pa or £^^''P'^ as well as the corre- 
sponding spacelike alternating tensor e'^"'' are all physically well defined in the sense of 
being independent of the choice of Aristotelian frame, since the effect on them of the 
boost transformation specified by (12) will be given trivially by 

t, = t, = Y'' , (17) 

and 

—pupa —pvpa 1 '-- ) ^ c . V 

Within the full group constituted by the (in general non linear) Milne transfor- 
mations characterised by (12) and (13), there is a linear subclass constituting the well 
known Galilei group, that is characterised by the requirement of preservation of the con- 
nection and the associated covariant differentiation operator V, for which the necessary 
and sufficient condition is that the boost acceleration vector should vanish. However 
for a generic Milne transformation the covariant differentiation operator will undergo a 
non trivial transformation, V i-^ V, specified by a corresponding transformation of the 
connection that will be given, using the definition (14), by the formula 

T;,^f;^ = T;^-t,a%, , (19) 

which has the noteworthy property of preserving the trace of the connection, i.e. we 
shall have 

r/p = r/, , (20) 

with the implication that for the evaluation of the simple divergence of a contravariant 
vector field, such as the displacement field introduced below, it will not matter which 
connection we use, i.e. we shall have V^^^ = V^^'^ 

Instead of working with the kind of fiat but ether gauge dependent connection that 
has just been described, it is useful for some purposes to work instead with a curved 
but gauge independent connection of a gravitational field dependent kind, that was first 
introduced by Cartan, and that is described in the preceding work referred to above [7]. 
However such a Cartan connection will not be needed in the present work. 

3. Relativistic correspondence 

A Newtonian space time structure of the kind described in the preceding section 
can be obtained as a low velocity limit from a corresponding relativistic theory on 
the supposition that the latter is approximately describable, in terms of an adjustable 
parameter c, by a Lorentz signature metric dS^ = g dx^ dx'^ , having the form 

gpu^Vnu-c'^tf.U, (21) 
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while, according to the preceding relations (5), the corresponding contravariant metric 
will be given by 

1 e^gv. (22) 

The quantity c in these expressions is interpretable as representing the speed of light 
with respect to coordinates of the standard Aristotelian kind as described in the pre- 
ceding section. In the absence of gravitational effects the metric (21) can be taken to 
be of flat Minkowski type, as given by a fixed value for c, but to allow for the effect of a 
Newtonian gravitational potential, (p say, it is necessary to take it to be given in terms 
of a fixed asymptotic value, c, by the formula 

= + 20 . (23) 

The degenerate Newtonian structure of the preceding section is then obtained by taking 
the limit c — > oo, which evidently gives 

^/^-^^M-, c-%.^-V.. (24) 

It is to be remarked that although the spacetime metric itself is degenerate in this 
Newtonian limit, the associated Riemannian connection, 

f = r''{9a(,,u) - \9,.,a) , (25) 

will be well behaved, and that it will agree with the fiat connection F^'^j^ of the preceding 
section in the absence of a gravitational field, i.e. when the potential (p in (23) is uniform. 
There will however be a difference in the generic case, for which it can be seen that the 
large c limit will be given by the relation 

f ^ r/, + t^t.7">,. ■ (26) 

This shows that in the Newtonian hmit F^''^ goes over to the Newton-Cartan con- 
nection that was denoted by cu^Pj^ in the preceding work [7]. This means that the 
associated Riemannian covariant differentiation operator will go over, not to the 
usual flat space differentiation operator Vi, of the preceding section, but to the Cartan 
type differentiation operator that was denoted in the preceding work [7] by L*^. 

In Newtonian dynamical theory, the concept of mass conservation plays an essential 
role. In relativistic theory, mass as such is not in general conserved, but in relevant 
applications it will nevertheless be possible to attribute most of the mass to other 
effectively conserved currents (e.g. that of baryons in a typical astrophysical context, 
or those of separate chemical elements in a typical non-nuclear terrestrial context) . Such 
conserved currents can be endowed with suitable mass weighting factors - e.g. the rest 
mass, m say, of an isolated proton or of a neutral hydrogen atom in the baryonic case 
- so as to provide what is needed in a Newtonian limit. A conserved current can 
be represented - without reference to any spacetime metric structure - as a Cartan 
type 3-form, with antisymmetric components n*^^ say, that is closed in the sense of 
having vanishing exterior derivative. This closure condition will be equivalent to that 
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of vanishing of the divergence of the corresponding current 4- vector that is given by the 
(Hodge type duahty) ansatz 

- j.e^'''"' , (27) 

in terms of the antisymmetric measure tensor associated with the spacetime metric g^j_y^ 
as given in terms of the modulus of the metric determinant \g\ by the condition that its 
components should be given by +||^||^'^^ or — whenever the indices are respectively 
even or odd permutations of the sequence {0, 1, 2, 3}. This means that it will be related 
to the non-relativistic spacetime measure of the preceding section by e^upa — cs^^p^ 
and £^^''p'^ — cr^e'^^'"^ , so that in the hmit c ^ oo one obtains n*^^ ^ Lf_iupa^'^)- 

A multiconstituent system may involve several constituent currents n^, which need 
not be separately conserved, but that combine to give a locally conserved total 



= E^n^ , V^vT = . (28) 

For a confined system, the corresponding globally conserved mass integral, M say, 
associated with a spacelike hypersurface - as specified by a fixed value of some suitable 
time coordinate - will be given in terms of the other coordinates and of the relevant 
mass parameter m by 

M^mJ n*^3 dx' dx' dx' , 71*^3 = c-X,,,n° . (29) 

In the applications under consideration, the evolution equations of the relevant cur- 
rents will be obtainable from a relativistic action principle of the world-line variational 
kind that is indispensible [15] for treatment of a constituent that is solid, and that 
is very suitable [16, 1, 2] (though other - e.g. Clebsch type [17] - possibilities exist) 
for treatment of a medium in which the relevant constituents are fluid. This kind of 
variational principle is based on a relativistic action integral 

jAc-^e,^^odx'dx^dx''dx\ (30) 

for which the action density will be decomposible in the form 

A = A,, + A,„, , (31) 

in which Aj^^^ is a relatively small intrinsic contribution. The dominant extrinsic contri- 
bution A^^j is the ballistic part - which is all that would be needed for the case of force 
free geodesic motion - namely the negative of the sum of the relevant rest mass-energy 
contributions, as given by 

Abai = -"^c^S^^x, (32) 

where m is the appropriate mass weighting factor and, for each constituent, is the 

corresponding number density as evaluated in the relevant rest frame. The frame in 
question is characterised by the corresponding timelike frame vector that is specified, 
subject to the normalisation conditions u^u^i, = — c^, by the expressions 

1/2 



< = < , ^x = c ^ (-<^xi^) ■ (33) 
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To obtain the Newtonian limit in which we are ultimately interested, with respect 
to the gauge specified by some chosen ether vector e^, it is convenient to express each 
frame vector in terms of a corresponding purely spacelike 3- velocity vector in the 
form 

K^u^{e^ + v^), <i. = 0, (34) 

in which it can be seen that the required normalisation factor is identifiable as the 

the time component of the frame vector with respect to coordinates of the standard 
Aristotelian-Cartesian type — x, — y, — z, x° — t, for which the metric 
takes the familiar form 

ds" = dx^ + dy"" + dz^ - c'dt^ . (35) 

This particular time component will evidently be expressible with respect to an arbi- 
trary coordinate system by the covariant formula 

= < = c (c^ - v^) , vl = < v^, = r^,,v^v^^ . (36) 
For the current itself we obtain the corresponding expression 

K = i^" + <) ' = ^x^° = K , (37) 

in which the relation between the rest frame number density and the corresponding 
ether frame component n° will be given by 

n.-Kc-'{c'-viy^\ (38) 



so that we shall have 

<-K' = <{^-"f)- (39) 
We thereby obtain a decomposition of the form 

A,, = A... + A_ , (40) 
in which the extrinsic action contribution is given by 

Ae.t = rn ^^^^ , 41 

which is evidently well behaved in the Newtonian limit, while the remaining contribu- 
tion will be expressible in terms of the total current (28) in the simple form 

Km = -m c'^n° , n° = n^^ . (42) 
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Although this remainder A^^_^^ will evidently be divergent in the large c limit, it can 
be seen that this does not matter from the point of view of the variational principle 
since the corresponding integrated action contribution, as given in terms of the standard 
coordinates used for (35) by 

^.e. = / K^c-%,,,dxdydzdt, (43) 

will be expressible in terms of the global mass function (29) simply as 

I^^ = -[Mc'dt. (44) 



In the kind of application under consideration, the admissible variations must respect 
the relevant current conservation law (28), so that they will have no effect on the global 
mass function M, which will therefore just be a constant that may be taken outside 
the integration. This means that the remainder term (44) will have a fixed value, so 
that from the point of view of the variation principle its inclusion is entirely redundant. 
The ensuing theory is thus unaffected by replacing the original relativistic Lagrangian 
density A by a new but equivalent Lagrangian density 

A =A-A =A,+A. (45) 

new rem ext ' int ' \ / 

from which the divergent term has simply been subtracted out. The new version 
A^^^ has the technical disadvantage of being gauge dependent, since its specification 
depends on the choice of the ether frame characterised by the vector e'*, but it has the 
important advantage of remaining well behaved in the Newtonian limit. With the usual 
assumption that (with respect to the chosen ether frame) the relevant squared space 
velocities and the potential are of the same order of smallness compared with the speed 
of light, we shall obtain 

as c — > oo. It can thereby be seen that in this limit the new version of the action 
density will take the form 

A.ew = A + 0{^}, (47) 

in which A is of purely Newtonian type, taking the standard generic form that was 
postulated at the outset in the preceding work [7, 8], namely 

A = Aint + Ae., , A^^, = A,^„ + A^„, , (48) 

in which the kinetic and gravitational potential contributions have their usual Newto- 
nian form 

Ah„ = 1^ ^x^xt^x ' Kot = -0 m E^n^ , (49) 
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4. Noetherian construction of stress - energy tensor. 



For a variationally formulated theory in a general relativistic context, the corre- 
sponding stress momentum energy tensor is commonly constructed directly by differen- 
tiation of the relevant action density with respect to the spacetime metric. The purpose 
of the present section is to describe the adaptation of such a procedure to the technically 
more complicated Newtonian case. 

In the preceding work on fluid systems [16, 7, 8] it was convenient to formulate the 
action in terms of physical fields (such as current 4-vectors) that are not entirely free 
but partially restrained as far as the application of the variation principle is concerned. 
However in the case of the solid systems with which the present article will be concerned 
it will be more convenient (albeit at the expense of greater gauge dependence) to work 
just with fields whose variation is unrestrained. The advantage of using an unrestrained 
variational formulation as we shall do here is that for an unperturbed field configura- 
tion satisfying the dynamical equations provided by the variational principle, the most 
general variation of the action will be unaffected, modulo the addition of a variationally 
irrelevant divergence, by the variations of the relevant dynamical fields, and will there- 
fore be given, modulo such a divergence, just by the contributions from variations of the 
given background fields characterising the relevant Newtonian or relativistic spacetime 
structure. 

In the relativistic case, the only independently given background field needed for 
this purpose is the metric gf^u itself. Provided that the other dynamical fields in the 
Lagrangian A obey the corresponding variational field equations, the generic action 
variation will be given by an expression of the simple form 



using the symbol = to denote equivalence modulo a divergence, i.e. modulo a term that 
is variationally irrelevant because its integral for a perturbation in a confined spacetime 
domain will vanish by Green's theorem. 

As has long been well known in the context of general relativistic theory [14], and 
as has more recently been demonstrated in the Newtonian case [8], the use of a fully 
covariant formulation makes it possible to derive useful Noether type identities by 
considering variations of the trivial kind generated by an arbitrary displacement field, 
say. This means that the variation of each (background or dynamical) field variable 
will be given by the negative of its Lie derivative. In the relativistic case, the relevant 
variations will be given simply by 



5A 



dA 



hi 



(50) 



'.V 



-SA = f£ A = r V^A ^ -AV^r , 



(51) 




(52) 



(53) 
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in which the relevant stress momentum energy density tensor can be read out as 



Ti"^ ^ 2-^ + Kg"'' . (54) 

By again removing a variationally irrelevant divergence, (53) can be rewritten equiva- 
lently as 

e^.f\ = , (55) 

which, since it must hold for a locally arbitrary vector field , shows that the variational 
field equations entail - as a generic consequence - a divergence condition of the well 
known form 

%f\^Q, (56) 

which in the absence of gravity, i.e. in a fiat Minkowski background, is interpretable 
as an energy-momentum conservation law. 

The (automatically symmetric) "geometric" kind of stress-energy tensor (54) needs 
to be distinguished from the (related, but in general different and not necessarily sym- 
metric) kind of stress-energy tensor commonly referred to as "canonical", which is 
constructed by differentiation not with respect to the metric, but with respect to the 
other dynamically relevant fields. Even in a special relativistic context, i.e. when grav- 
itation is negligible so that the physical metric can be taken to be fixed, it is still 
perfectly legitimate (though the possibility of doing so is not widely realised) to exploit 
the greater convenience of the "geometric" construction via the consideration of virtual 
( "off shell" ) variations of the metric. However the construction of a "geometric" stress 
energy tensor is not quite so straightforward in a non-relativistic Newtonian framework, 
due to the degeneracy of the metric, which makes it harder to avoid the inelegancies of 
the traditional "canonical" treatment. 

Although it is not quite so simple and convenient as in the relativistic case, a "ge- 
ometric" type ansatz for the construction of an appropriate variational stress-energy 
tensor can nevertheless be obtained in a Newtonian framework using the 4-dimensionally 
covariant formalism set up in the previous section. Such an ansatz was developed in the 
preceding work [8] , where it was shown how, in the case of a simple or multiconstituent 
fiuid model, the action density function A for a 4-dimensionally covariant variational 
formulation provides a Noether identity that leads automatically to a natural "geomet- 
ric" type ansatz for a corrcponding non-relativistic strcss-momentum-energy density 
tensor T^^. The present section describes the way to obtain the appropriate Noethe- 
rian ansatz for the non-relativistic "geometric" stress-energy tensor in a manner that 
is rather simpler than was required for the partially restrained case dealt with in the 
preceeding work [8]. 

The kind of system under consideration is one governed by a non relativistic action 
integral of the form 

1= j Ae^^^^dx^dx^dx^'dx^, (57) 

in which the action density A is a scalar of the generic form (65) that is formulated as 
a function just of the relevant (variationally unrestrained) dynamical fields and their 
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gradients, and of the background fields and characterising the Milne structure of 
Newtonian spacctimc, as well as on a gauge dependent ether frame vector field e'^. The 
involvement of the latter will violate the Milne and even Galileian invariance of the local 
action density, but will not violate the required covariance of the ensuing field equations 
provided the effect of an ether gauge transformation A i— > A is just to add on a pure 
divergence term, which will not contribute to the integral of a local variation. This 
requirement is conveniently expressible as A = A, again using the symbol = to denote 
equivalence modulo a (variationally irrelevant) divergence, for whose construction (due 
to the ether frame invariance of the measure tensor e^i,p^ as remarked at the end of 
the preceding section) it makes no difference whether we use the original covariant 
differentiation operator, V, or the modified operator V. 

As remarked above, the simplification in the unrestrained case dealt with here is 
that for an unperturbed field configuration satisfying the dynamical equations provided 
by the variational principle, the most general variation of the action will be unaffected, 
modulo the addition of a variationally irrelevant divergence, by the variations of the 
relevant dynamical fields. Modulo such a divergence, the local action variation will 
therefore be given, in the non-relativistic case, just by the contributions from variations 
of the uniform background fields 7^*^, and e^^, as well as of a generically non-uniform 
gravitational potential (p in cases for which the latter is introduced as given non dynam- 
ical background. The variation will therefore be given, modulo a divergence, in terms 
just of a set of tensorial coefficients dh./d'^^" , dA/dtf^, and dA/de'^ (of which the latter 
would not be needed in a model with ether gauge independent action density) by an 
expression of the form 

dA OA aA dA 



It is to be remarked that this expression is not by itself sufficient to fully determine the 
coefficients involved. Further suitably chosen conventions, of which the most obviously 
appropriate is that the tensor dA/d'-f'^'^ should be symmetric, in view of the algebraic 
restrictions on the independence of the variations involved, which by (1) and (5) must 
must evidently satisfy 

t^^Sj'^" = -r^St, , t,5e'' = -e'^St^ . (59) 

For the actual evaluation of the coefficients in (58) it will also be useful to have formulae 
for the variations of other associated spacetime background fields that may be involved 
in the explicit formulation of the action, such as the covariant spacetime metric whose 
variation will be given by the formula 

^Vf^i^ = -Vf^pVuah'"' - 2t(^r]^)xSe^ , (60) 
and the spacetime alternating tensor whose variation will be given by the formula 

(in which the ether frame dependence of the two separate terms can be seen to cancel 
out to give a gauge invariant total) . 
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As in the relativistic case above, we now consider the effect of variations of the 
trivial kind generated by an arbitrary displacement field, say, so that the variation 
of each (background or dynamical) field variable will be given by the negative of its Lie 
derivative. The relevant formulae are thus given by 

-5A^i£A = e^,A, (62) 

-S-f^"" = C£ 7'"" = e^pl'"" - 27" ^"VpC^ , (63) 

-St^ = ^£ = eVpt, + tpV^e , (64) 

-Se" = i£ e'' = i^Vpe" - e^WpC" ■ (65) 

-5<f)^i£,<P = e^p<P, (66) 
The first of these formulae can be rewritten, modulo a divergence, as 

SA^AVpe, (67) 

while in view of the uniformity properties (4) and (6) of the unperturbed background 
fields, the next three will reduce to the form 

5^^''' ^2Y'^^VpC\ (68) 

5ti,^-tpV^e, (69) 
Se^" = e'Vp^ . (70) 
For such a displacement variation, the relation (58) will therefore reduce to the form 

nv;.r = prv.0, (71) 

in terms of a stress-momentum energy density tensor T'^ a gravitational mass density 
p that can be read out as 

n = A^^-2|^r''+|^t.-|^.". (72) 

and 

By a further divergence adjustment the equivalence (71) can be rewritten as 

r(v^n + pv.0)-o, (74) 

which means that for a displacement confined to a localised spacetime region outside 
which the hypersurface contribution provided via Green's theorem by the unspecified 
divergence term will vanish) we shall have 



/ 



r(V^r't + pV,0) £0123 dx' (75) 
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Since this identity must hold for an arbitrarily chosen displacement field in the 
spacetime region under consideration, it follows that we must have 

V^n = -pV.(/., (76) 

as a Noether type identity. We have thus established a theorem to the effect that 
the conservation law (76) will hold automatically for the geometric energy momentum 
tensor obtained from the covariantly formulated Newtonian action density A accord- 
ing to the prescription (72) whenever the dynamical field equations provided by the 
corresponding unrestrained variation principle are satisfied. 

In the simple applications to be considered below in the present article, the local 
energy momentum conservation law (76) will, by itself amount to complete set of dy- 
namical field equations, but it will only be a subset theroff in more general cases, such 
as the multiconstituent applications that we plan to deal with in subsequent work. 

5. The Material Projection concept 

The historical development of the standard textbook theory (see e.g. Landau and 
Lifshitz [18]) of a perfectly elastic solid in a Newtonian context is attributable to many 
people, among whom Cauchy is perhaps the most notable. However, as in the multicon- 
stituent fiuid case [5, 7], the insight needed for the formulation of a generally covariant 
version of the theory comes rather from its relativistic generalisation, for which a fully 
covariant formulation has always been used, as an indispensible necessity, from the 
outset. Some of the earliest work on the appropriate relativistic theory of a simple 
perfect solid was carried out in a purely mathematical context by Souriau [19], and 
by DeWitt [20] (who needed it as a toy model for testing techniques to be used in 
the quantisation of gravity). Its development as a realistic physical theory for use in 
the kind of astrophysical context (particular that of neutron stars) that motivates the 
present work was initiated (in the aftermath of the discovery of pulsars) with Quintana 
by one of the present authors [21], in a formulation [22] that was shown to be elegantly 
obtainable by a variational approach [15, 16] that will be used as a guide for the present 
work, whose purpose is to derive its Newtonian analogue. 

The material projection is based on the simple consideration that the intrinsic struc- 
ture of a solid is essentilly 3 dimensional. This means that the mathematical entities 
(such as differential forms) that will be used in a variational principle governing the 
dynamic behaviour of the solid should be defined over a 3 dimensional manifold, X 
say. The prototypical example of such field is the elastic-stress tensor, whose definition 
should not depend on the solid's history, whereas its explicit value obviously will do so. 
The requisite 3 manifold X is identifiable as the quotient of spacetime by the worldines 
of the idealised particles (representing microscopic lattice sites) constituing the solid , 
so that each point on X can be considered as the projection of the world-line describing 
the extrinsic motion of the relevant particle. Thus a patch of local coordinates, let us 
say (for A=l,2,3), on X will induce a corresponding set of scalar fields g-^ that will 
be given as functions of the local spacetime coordinates, (for = 0, 1, 2, 3), on 

The local scalar fields can be interpreted as a set of comoving - Lagrange type - 
coordinates on spacetime. They might even be used to specify the choice of the space 
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coordinates by taking — q^, x'^ = x^ = q^. However such a choice is more hkely 
to be convenient in a general relativistic context than in that of a flat Minkowski back- 
ground, or in the Newtonian case, for which a more commonly convenient choice is to 
use background coordinates x^ fiat (respectively Minkowski or Aristotelian-Cartesian) 
in order to simplify the procedure of covariant differentiation (in a manner that is 
not possible in the general relativistic case) by setting the connection coefficients F^^^ 
everywhere to zero. 

Since the whole worldline of the particle is mapped into a point in X the g"^, when 
viewed as scalar fields on A4, will obviously be characterised by the property 

uLq"" = w'^V^g^ = w'^gf^ = , (77) 

where the tangent vector field to the worldline is subject to the standard normali- 
sation condition given by u^u^ = —(? in the relativistic case, and hence by = 1 in 
the Newtonian limit (24). Using the symmetry property (2) of the connection, which 
ensures that we shall have V^g^^ = Vjygj^, the relation (77) in turns implies that 

u£,q^^ = u^V^q^^ + gf, VX = . (78) 



Let us now consider the example of a material 1-form on the manifold A", with 
components say. When pulled back via the material projection, this material form 
induces a spacetime covariant 1-form on the spacetime manifold Ai with components 
An given simply by 

^m = ^bC' (79) 
with the direct implication that one will have 

A^u^" = , u£,An = IvTVi^A^^ = (80) 

for any such material 1-form. Conversely if any spacetime 1-form is such that it satisfies 
both of the conditions above, then it is uniquely determined by a material 1-form 
through the puUback operation. 

The generic defining property of the kind of simple perfectly elastic model to be 
considered here is that the action should depend only on the rheometric position coor- 
dinates g^ and on the corresponding induced contravariant metric components 7^^, or 
equivalently on the corresponding contravariant components 7-*^ , which are defineable 
by the reciprocity relation 

r^iBc = K ' (81) 

and which will be given in the relativistic case simply as the rheometric projection of 
the contravariant spacetime metric, 

7-^"=5'^'gX-, (82) 
while in the Newtonian limit (24) the required components will be given by 

7""=7'^'^9fX-- (83) 
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It is to be remarked that, in solid models of the kind appropriate for typical (low 
pressure) laboratory type terrestrial applications, the rest-frame energy per particle will 
commonly have an absolute minimum for some preferred value, k.^^ say, of the induced 
metric ■j^g. In such such cases K^g can be considered, and exploited, as a natural 
fixed (positive definite) Riemannian metric on the material base space. In other words 
K^gdq^dq'^ represents the "would be" relaxed distance between the chosen particle and 
a nearby one in the frame of the former, assuming the attainability of such relaxed 
state (which could be defined by extracting a piece of the continuous medium - i.e. 
a neighbourhood of the considered particle - and letting it reach the minimum local 
energy density state). In generic circumstances however, such a preferred rheometric 
metric may be ill defined, since a local state of minimised energy need not exist. This 
caveat applies particular in cases of very high pressure (such as will occur in deep stellar 
interiors) from which a process relaxation might lead, not to a minimised energy state 
for the solid, but merely to its vaporisation as a gas. 

The formula (83), while defined in the Newtonian limit, will also be vahd in the 
relativistic case provided it is interpreted in terms of the relevant rank-3 worldline- 
tangential projection tensor, which will be given by 

r^^S^ + c-^^u,, (84) 

in the relativistic case, and which will will go over, in the Newtonian limit (24), to 

7^ = 5^:-«%. (85) 

(Except in the case of a static configuration - for which the ether vector e'^ may be 
chosen to coincide with -u^ - the variable fiow-tangential projection tensor defined by 
(85) must be distinguished from the uniform ether frame dependent projection tensor 
(5) that was denoted by the same symbol in the preceding work [7, 8] but that is denoted 
here by 77^). 

We shall use the induced base metric of (81) for raising and lowering of the material 
base indices in the usual way, as illustrated for a covector with components by 

A- = r'^Ag , A, = 7,,A- , (86) 



The (variable) covariant metric as defined by 

7m<^ = TAB^fM^fr.- (87) 

can be used in conjunction with the (uniform) contravariant metric 7^^^ for an unam- 
biguously reversible index raising lowering and raising operation for vectors and tensors 
that are orthogonal to the material fiow in the manner exemplified, according to (80), 
by the pull back (79) of A^ namely 

A^^A^q^^^-f^,A\ (88) 

whose raised version will project onto A^: 

A^'^Y'^A^, AH^^O, q^'A^'^A^. (89) 
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Is to be noted that although the procedure defined by (86) is invariant with respect 
to the choice of ether frame, it does depend on time. Thus if is taken to be a fixed 
vector potential characterising the "frozen-in" magnetic field that will be introduced 
below, so that its time derivative will vanish, the corresponding contravariant vector 
will neverthess be time dependent: 

A^ = Q ^ = Y'^^B ■ (90) 

It can be seen from (4) and (78) that the time derivative, along the worldhnes, of 
the induced base metric will be given by 

Iab = ^Oab = -^acIbo^'' , = -2^-- = -Sgf^gf.^'^'^ , (91) 

in terms of the (symmetric) strain rate tensor 

Qfxv ^ Y^'^VpU''^ . (92) 

which will automatically satisfy the orthogonality condition d^^u^^ = in the relativistic 
case, so that we shall have O^^t^ = in the Newtonian limit. 

The time derivation considered so far has concerned only quantities such as the base 
space components that have the status of scalars from the point of view of the back- 
ground space time. We now extend the dot notation to quantities that are tensorial with 
respect to the backkground spacetime by defining it to indicate covarant differentiation 
with respect to time along the world lines, meaning that ' will indicate the effect of 
the operator u^Vp in the manner illustrated by the definition of the acceleration vector 
which will be given by 

uf" = u^VpU^" . (93) 

so that it will satisfy the condition ii'^u^ = in the relativistic case and iiH^ = in 
the Newtonian limit. It can be seen that the corresponding covariant time derivatives 
7^,^ = M''Vp7^jy of the gauge invariant metric fields 7^,^ defined by (92) will be given in 
the relativistic case by 

2 

inv = -^u^fiUu) , (94) 

and in the Newtonian limit by 

= -2t^f^7u)pu' , Y = -t.u'' . (95) 



It is useful for many purposes - and will be indispensible for the discussion of the 
Newtonian limit - to introduce an appropriate fixed measure form, with antisymmetric 
components N^^^ — -/V,^^^] say, on the rheometric base manifold. Such a measure 
will typically be interpretable as representing the density of microscopic lattice points 
in an underlying cristal structure. Such a measure will determine a corresponding 
contravariant spacetime current of the kind introduced in (27) by a pull back relation 
of the form 

^.Bc9f/.9" (96) 



18 



The corresponding scalar number density n in the rest frame of the medium will evi- 
dently be given by 

= I^NascN^^f r^r^r' ■ (97) 

It is then obvious that the current density so defined will be automatically conserved 
due to the closure of the 3 form N^^^ (being a 3 form on a 3 manifold) : 



V[AiV^,p]=0, VX = 0. (98) 

In the simple purely elastic application considered here, the formalism defined above 
will be enough to describe the dynamics of the solid. However such a simplification will 
not be possible in the more general applications envisaged for future work, such as 
transfusive exchange [2] of matter between distinct chemical constituents, in which one 
no longer has conservation of the solid's current. 

6. Action for a simple perfectly elastic solid 

Whereas the internal energy depends only on the density when a simple fiuid is 
considered, in the case of a solid it will also be necessary to allow for the reaction of 
the internal energy not just to changes of volume, but also to changes in shear strain. 
This is done by allowing Aint to depend on which, because of its time dependence 
can be interpreted as a Cauchy type strain tensor. This kind of (perfectly elastic) solid 
model includes the category of simple (barotropic) perfect fiuid models as the special 
case for which the dependence is only on the determinant of j^^. For a medium that is 
perfectly elastic in this sense, the generic action variation will be given in the relativistic 
case by the formula 

c-T , (9A , 

SA = Sj^'^ + T^Sq^ , 99 

which characterises the partial derivative components needed for the specification of 
the rheometric stress tensor S^^g and its spacetime puUback 

dA 

S,. = S,,q-^q-,, S,,=2-^-A^,,. (100) 

It can be seen from (54) and (84) that the complete stress energy tensor will be given 
in terms of this quantity by an expression of the standard form 

f"" = puV - S"" , (101) 

in which comoving mass-density density p is given by the simple proportionality relation 

A=-pc^. (102) 

This is all that is needed for the formulation of the variational field equations, which 
in this case are given completely just by the Noetherian divergence condition (56). 

In typical physical applications the mass-energy function in (102) will have a mini- 
mum value Pj,^^ - corresponding to a maximum value A^^ of the elastic action function 
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A - for any given value of the determinant I7I of the induced metric 7^^, or equivalently 
for any given value of the number density n as specified by (97). This means that the 
action will be decomposible in the form 

A = Aflu + Ais' (103) 

in which the - usually rather small - remainder term - A^j^ is a negative indefinite rigid- 
ity contribution contribution without which the medium would be of purely fluid type. 
The - usually dominant - perfect fluid contribution A^,^ will itself be itself decomposible 
in the form 

= + A,.. , (104) 

in which Ap^.^ is a pressure energy contribution that in typical applications will also 
be relatively small compared with a dominant ballistic contribution Aj^^j given by an 
appropriate choice of the mass parameter m in the general formula (32) which, in the 
single constituent case considered here, reduce to the trivial form 

Abal = -PC^ p^mn, (105) 

so that the mass density (102) will be expressible in the form 

p^p + c-^€, (106) 

with the relatively small internal energy contribution given according to (31) by 

^ = -A.. • (107) 

The action decomposition (103) corresponds to an energy density decomposition of the 
form 

6' = g' +6*. , (108) 

pre ' rig ' V / 

in which the rigidity contribution will vanish when & is minimised for a given 
value of n. The non trivial pressure term 

= -^P.e , (109) 

in (104) will evidently be a function just of the undifferentiated base space coordinates 
and of the scalar density n, which means that, as the analogue of (99), its generic 
variation will be given by 

with 

Tl 

= ■ (111) 

It is evident from (100) that the ballistic term will make no contribution at all to the 
strain tensor, while the contribution from the pressure energy term will of course be 
purely isotropic: 

QHV^r^ Q iiv ^ _p tiv p P£l _ ^ fllO'l 
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This means that the purely intrinsic action contribution given, according to (108), by 

A. = A + A. , (113) 

int pre rig ' \ / 

will be sufficient by itself to determine the the complete stress tensor in (101), which 
will take the form 

int pre ' rig ^ ' 

It is important to notice that in this simple elastic case the effect of a worldline 
displacement along the worldlines will have no effect so substitution of for in (55) 
will merely give an identity 

li'^V.T^^O. (115) 

A complete set of dynamical equations will therefore be provided just by the orthogonal 
projection of (56) which will be expressible as 

(p7t - c-'S^,) ii^ = i^r^'^SP, . (116) 



7. Newtonian dynamics of a simple perfectly elastic solid. 

To obtain the Newtonian limit of the simple relativistic elasticity model charac- 
terised by the stress energy tensor (101) that was set up in the preceeding section, it is 
now straightforward to apply the procedure described in Section 3, according to which 
we should use a Newtonian action integral 

1= j Ive^^^^dx^dx^dx^dx\ (117) 

with a Lagrangian of the form (48), which in this simple elastic case reduces just to 

A = A,.„ + Ap..-€', (118) 

with an internal energy function € of the same form as in the relativistic case, and 
with the Newtonian kinetic and potential energy terms given by expressions of the same 
form as for a simple fluid, namely 

K. = \P^\ Apot = -#, (119) 
in which the 3-velocity v'' is charcterised in terms of the chosen ether frame frame, 

by 

v^^ = u^'-e^' , = r]^^ u"v" , (120) 

and relevant Newtonian mass density will be given simply by 

p — mn , (121) 
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for a number density n that, in accordance with (96), will be expressible in this limit 
as 

n = nX = l^e^^'n^., , = nu^ = ^£'^^'"^<p. ■ (122) 

The current 3 form n*^^. itself is obtained from the fixed 3- form A^^^^ on the material 
base space by the construction (96), which automatically ensures the satisfaction of the 
non-relativistic conservation law in its usual form 

VX = 0. (123) 

It is to be remarked that, as in the fluid case [7], the gravitational coupling term Ap^^ 
will be unaffected by linear (Galilean) gauge transformations but not by accelerated 
Milne transformations, while the kinetic term A^.^ is of course even more highly frame 
dependent (not even Galilei invariant). For the purpose of comparison with the preceding 
work, the terms in the Newtonian action can usefully be regrouped in the standard form 

"^mat "^pot "^mat "^kin ^ ' ("^^^) 



In the perfect fluid case this elastic energy density & will be given on each worldline 
(as specified by the values q^) as a, corresponding function (in the barotropic case [7] 
everywhere the same function) of the number density n, which can be seen from (121) 
to be given as a function of the scalar product fields 7^^ by the determinant formula 

= l^N.^cN^^.r^r^r^ ■ (125) 

The generalisation from a perfect fluid to a perfectly elastic solid is made simply by 
taking to be a generic (worldhne dependent) function of the scalar product fields 7^^ 
(i.e. of the tensorial projection of 7^^^ onto X) instead of being restricted to the depend 
just on the determinental combination (125) as in the fluid case. This means allowing 
to be affected not just by changes of volume but also by shearing strain, whose effect 
is allowed for by supplementing the purely fluid contribution g'pre in (107) by the extra 
rigidity term €^rig- 

Using (60) and (61) we obtain the formulae 

Su^" = -u^'u^'SU (126) 

and 

for variations in which the fields are held constant - meaning that the world lines 
remain fixed - it can be seen from the defining ansatz (72) that the kinetic and gravi- 
tational potential contributions to the stress momentum energy density tensor will be 
given by expressions that are already familiar from experience [7] with the simple fluid 
case, namely 
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in which the frame dependent 4- momentum per particle is given, using the notation 
of (120) by 

P^^ = m{v^ - '-vX) (129) 

and 

TJ. - -<l>p''t. , ff = mn^. (130) 

The corresponding stress tensor S*-^^ is definable in the traditional manner in terms 
of the effect on the energy per particle €/n of an infinitesimal strain variation 5'-)^^ 
according to a prescription of the form 

(131) 

By comparing this with the formula (111) for 5n, it can be seen that this contravariant 
stress tensor will be expressible independently of the number density as 

3f 

^ 2- + 6'7^« , (132) 

which means that its covariant (index lowered) version will be given by 

^.. = -2^ + 6^7... (133) 

In the manner described above this base tensor will determine a corresponding space 
stress tensor 

S^.u = S^.q^'^^q'',, , S^^u" = , (134) 

in terms of which the rate of change of the energy density along the world lines will be 
given by 

e = S^,0^- - €9 , (135) 
where 9 is the expansion rate, as given by 

9 = Tm-^''" = V^ii^ = -fi/n . (136) 

Using the formula 

g^AB 

= (137) 
obtained from (83), it can be now seen from (131) that we shall have 

^ = \{€;^AB - S,,)q^X^ = i(€^7.. - V) • (138) 
The relation (107) thus gives 

dA. . 



US^.-€^^,), (139) 
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so that by the general formula (72) (since no dependence on e'^ or is involved in this 
case) the internal contribution to the stress momentum energy tensor will be given by 

= + - 5p.) , (140) 

which simplifies, using (85) and the (orthogonal) index raising operation exemplified 
by (88) to 

TJ, = -€un, - ■ (141) 

Thus, combining this gauge independent internal contribution with the ether frame 
dependent kinetic contribution (128), we end up with the complete material stress 
energy tensor 

T = T + T t , (142) 

in the form 

TZ. = ^'^^ - S% . (143) 

in which the relevant (gauge dependent) momentum per particle covector is given 
by 

m,=p,- {€/n)U = m[v,- {\v^ + € / p)U^ . (144) 

Since it is not orthogonal to the flow this momentum 1-form is not completely deter- 
mined just by the corresponding - purely kinematic - contravariant momentum covec- 
tor, namely 

ml" = = mvi" . (145) 

In the special case of the perfect fluid limit we shall simply have S'^^ — —P'-f'^", 
where P is the ordinary scalar pressure. It is important to note, in this case, that the 
momentum covector rrii, introduced here will not be quite the same as the material 
momentum covector p,^ that was introduced in the preceding work [7], since rrii, is 
defined in terms of the integrated internal energy per particle, namely 6*/n, whereas fi^ 
was defined in terms of the differential internal energy per particle, namely the chemical 
potential, x— {€ + P)/n. This means that, in the (barotropic) perfect fluid hmit, the 
relevant material momentum 1-form will be given by /i^ = — [P/n)tfj,. 

The complete stress energy tensor, including allowance for the gravitational back- 
ground, can now be obtained as 

T^-TZ.^TJ^u. (146) 

which by (130) flnally gives 

Tt = ri>'{m, - m(f>U) - 5'^ . (147) 

for the tensor whose divergence will provide the required dynamical equations according 
to the Noether relation (76). 

For a compound system (as exemplified by the muticonstituent fiuid models studied 
in the preceding work [8]) the 4-indcpcndent components of the Noether relation (76) 
would not by themselves be sufficient to fully determine the dynamical evolution. How- 
ever for a simple medium such as we are considering here, for which the dynamics are 
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completely describablc just in terms of the motion of the world lines, as given by the 
evolution of the 3 independent scalar fields g-*, the evolution is actually overdetermined 
by the 4 components of (76), whose contraction with 4- velocity, namely 

u'^f.^O, /, = (V^T't + pV.(/))=0, (148) 

merely gives the kinematic identity (134), which must be satisfied as a mathematical ne- 
cessity, even under circumstances more general than those considered here, in which the 
space components of the force density acting on the medium might not be zero. The 
underlying reason for this identity is that whereas the dynamical equations were needed 
for the derivation of the Noetherian condition (74) for an arbitrary displacement vector 
field it would evidently hold as a mere identity for a displacement along the fiow 
fines, i.e. for oc u'^, in the simple elastic case. It is to be remarked that the logic could 
be reversed as, was done in the cited work on the relativistic case, which started [21] 
by postulating the analogue of (148) as a condition needed for consistency, and then 
went on [15] to derive the action formulation as a consequence. 

The upshot is that the complete system of dynamical evolution equations for a 
simple elastic solid model will be given just by the 3 independent components of the 
space projection of (76) namely 

7'^'^(%n + pV,0) = O. (149) 

Using the conservation law (123), it can be seen from the formula (147) that in terms 
of the frame dependent gravitational field vector, 

^ -^/^-V,(/. (150) 

this system will be expressible in the form 

piii'^-gn^V^S'^'', (151) 

in which it is to be observed that each side (though not the separate terms on the left) 
is satisfactorily invariant, not just under linear Galilean tranformations, but even under 
arbitrarily accelerated Milne transformations. 
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8. Derivation of the cheiracteristic equation. 

As in the relativistic case [25] , let us now seek the conditions governing the a covector 
say that is normal to a characteristic hypersurface across which the relevant field 
quantities n, u^, and S^^^ have discontinuous gradients, using the standard method 
of Hadarmard, which exploits the condition that the discontinuity of the gradient of 
a continuous scalar field must be proportional to the normal covector A^, so that in 
particular for the density we shall have 

M=hX^, (152) 

for some corresponding scalar discontinuity amplitude n. The associated unit propaga- 
tion vector i/^ characterised by 

r^i^.^f. = 1 , ly^.W^O, (153) 

and the propagation velocity, v say, relative to the local rest frame, of the discontinuity 
are specifiable by taking A^ to have the standard normalisation so that it takes the form 

A;u = i^^u + V c~^w^ , (154) 

in the relativistic case, and hence 

>^fi^J^n- vtu, , (155) 

in the Newtonian limit. In terms of the same discontinuity covector as in (152) the 
discontinuity of the gradient of will be given by an expression of analogous form, 

[VX]=«"A^, (156) 

in terms of a corresponding vectorial discontinuity amplitude , while similarly for 
the stress tensor we shall have 

[V^^'^l^^^^. (157) 



Since the evolution of n and 5'''^ is kinematically determined by that of the flow 
lines, the corresponding gradient discontinuity amplitudes fi and S^^ will be determined 
by the velocity gradient discontinuity amplitude u. In the case of the number density 
n it can be seen from (136) that we shall have 

ui'V^n = n = -pO = -nV^u>' (158) 

so by taking the discontinuity we obtain 

u^'X^n = -nX^u^ . (159) 

The normalisation conditions u^u^ = in the relativistic case and u^t^j, = 1 in the 
Newtonian case imply corresponding restrictions u^u^ — 0, uH^ — respectively, with 
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the implication that for ?}/', as for z/^ we can unambigously and reversibly raise and 
lower the indices by contraction with 7''^ and 7^,/. It follows that (159) will reduce to 
the simple form 

vn = nv^u'^. (160) 

To write the corresponding relation for S^^ we need the relevant elasticity tensor, which 
is defined in such a way as to have the symmetry properties 



^ABCD — £JCDAB _ £](AB)(CD) ^ (161) 

by the ansatz 

JT'ABCD ^ Ar, { — 



^ABCD^^^ (_]=2- +,^^^7^^. (162) 

dlAsOlcD^n^ dlcD 



In terms of this highly symmetric elasticity tensor, the less highly symmetric Hadamard 
elasticity tensor that will be needed below is specifiable as 



\ABCD 



^ J^CDAB ^ ^ABCD ^ ^AC gBD _ (^gS) 



It follows from (162) that the time derivative of the stress tensor in the material 
base space will be given in terms of that of the strain tensor by 

S^B ^ ^(^^ABCD _ gAB^CD^A^^^ _ (164) 

For the purpose of evaluating the time derivatives of the corresponding space time 
tensors, contravariant base tensors are less convenient than the corresponding covari- 
ant tensors, whose time derivative can be seen, by (78), to puUback directly onto the 
corresponding Lie derivative in the manner illustrated in the case of the stress as 

?f^?" -^^B = ^£,S^, = ^''V,^^, + 2Sp^^V,)uP . (165) 

We therefore need the formula obtained by swapping covariant with contravariant in- 
dices in (164) which gives 

Sab = -'^i^ABcn " '^.bTc. + 4^c(.7., Jt'^" , (166) 

from which, by (91) we obtain 

= {E^,,, - S^^j,, + 45,(^7.),)^'"^ . (167) 

Combining this with (165) and using the definition (92) of the expansion rate tensor 
9'^'^ we obtain an evolution equation for the stress tensor in the form 

H^VpS-^. = -25,(^V,)«'' + (£;^./ - 5^.7/ + 2^p(/.7.^ + 25'^(^7,),) V.m" . (168) 

Taking the discontinuity of the gradients in this relation we obtain, as the analogue of 
(159), 

u'XpS^, = -2Sp^^X,)uP + (£;^,/ - 5^,7/ + 25,(^7,^ + 25'^(^7,),)\ii'' . (169) 
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After projecting out the time component by contraction with 7^^^. this leaves, as the 
analogue of (160), an expression giving the stress gradient discontinuity amplitude as 
a function of the velocity gradient discontinuity in the form 

V r'Y'Sp, = -{Ei^'^P, - S'^'^Ya + 2-S,('^7'^0^"^P, ■ (170) 

We now have all that is needed for processing the gradient discontinuity relation 
provided by the equations of motion, from which one obtains the dynamical equation 

pu''X,y^X-KY'r''Spa: (171) 
in terms of a tensor y'*^ that will be given in the relativistic case (116) by the relation 

pcYu = ipc' + €)Y, - 5t , (172) 

but the reduces in the Newtonian limit (151) simply to = 7^*^,. 

By substitution of the kinematic formula (170) into the dynamical condition (171) 
we finally obtain the required characteristic equation in the form 

(i;Vy'^'^-Q'^>. = 0. (173) 

This is an effectively 3-dimensional eigenvector equation with -u^, the square of the rela- 
tive propagation speed, as eigenvalue, for which the eigenvector is the covariant velocity 
gradient discontinuity amplitude as characterised by the orthogonality condition 

u^u^^Q. (174) 

In terms of the Hadamard elasticity tensor specified according to (163), the characteristc 
matrix Q^^ can be seen to be expressible as a function of the propagation direction, as 
indicated by the spacelike unit covector by the formula 

= A'^^'^VpZ/^ . (175) 

The simplest application of this formula is of course to the case of a medium that is 
intrinsically isotropic (as will typically be the case in macroscopic applications involving 
averaging over a large number of randomly oriented mesoscopic crystals) and that is in 
an undeformed, though perhaps highly compressed state, with energy density €q say. 
In such an undeformed state the stress tensor S^'^ will reduce to an undeformed value, 
Sq^ say, that will be characterised just by a pressure scalar Pq, in terms of which it will 
take the form 

S^"" = -PoY" ■ (176) 

It follows that the tensor y^'^ in (173) will reduce to an undeformed value y/^'^ given in 
the relativistic case by 

Vo''^ = (1 + 7- ■ (177) 
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As discussed in the cited work [21] on the relativistic case, for such a state the elasticity 
tensor E^^p" will reduce to a corresponding isotropic value Eq'^'^'^ that is expressible in 
the well known form 

E^"^ = (/3o - \Po)l^''l'^ + 2(/io + P^){l^^'r^'' - ll^^'in , (178) 

in which the coefficients (3q and /xq arc respectively interpretable as the bulk modulus 
and the modulus of rigidity In the particular case of a perfect fluid the rigidity will 
vanish, = 0, and the bulk modulus will be given by the derivative of the pressure with 
respect to fractional volume change, (3 — ndP/dn. According to (163) the Hadamard 
elasticity tensor A^^'^^"' reduce to a corresponding isotropic hmit value given by 

Ar"" = Po-r-r + 2p^y^^y^' + 2/xo(7''^''7"^" - Iy'Y") , (179) 

which is such that the antisymmetric pressure term will cancel out in the formula for 
the characteristic matrix Q^^ ^ leaving an expression of the same form Qq" as is familiar 
in the low pressure limit, namely 

Q^' = (/3o + \lJ^o>^i^' + IJ^oY' . (180) 
9. Faraday - Ampere magnetodynamics. 

So long as it acts merely as a given prescribed background, an electromagnetic field 
can be treated within a Newtonian framework in a manner that is satisfactorily Galilei 
and even Milne invariant [26]. However when it is necessary to treat it in its own 
right as an active dynamical field governed by an unrestricted electric current source 
then it is necessary to sacrifice Galilean (and hence a fortiori Milne) invariance, which 
was brutally violated by the introduction of a physically preferred ether in Maxwell's 
original formulation, but more elegantly replaced by Lorentz invariance in Einstein's 
special relativistic treatment. 

A satisfactorily Galilean and even Milne invariant formulation in terms of an anti- 
symmetric field 2-form with components F^jy = —Fy^ can however be set up as a self 
consistent approximation in cases for which only a subset of three "magnetic" degrees of 
freedom are dynamically independent, while the other three "electric" components (out 
of the six contained in F^,^) and all the components of the current are treated merely 
as passively derived fields in the manner that will be described immediately below in 
this section. 

In the following section it will be shown how such an ether gauge invariant Faraday- 
Ampere type model can be coupled, in a variational formulation, to a simple perfect 
barotropic fluid or elastic solid model, of the kind described in the previous section, 
in the special case of "perfect conductivity", meaning the case for which the field is 
"purely magnetic" in the sense that, with respect to the local rest frame specified by 
the velocity 4 vector of the fiuid, the relevant electric components are zero. We thus 
obtain a 4-dimensionally covariant formulation of the non - relativistic version of what 
is known in the fluid case [27] as a perfect magnetohydrodynamics. 
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In the generic case, the electric and magnetic fields and S^,^ can be defined, with 
respect to an ether vector e'^, by the decomposition 

F^u = + '2E[fj,t,y] , (181) 

subject to the conditions 

E^e^ = , B^^e" = , (182) 
which are equivalent to the specification 

E^^F^^e", (183) 

where, in the Newtonian case we are concerned with here, is the preferred time 
covector introduced in (1) (while in the relativistic case it would be given in terms of the 
spacetime metric by = —g^ivU^)- Under the action of an ether gauge transformation 
of the form (12), as generated by a spacelike boost vector field h^, these fields will 
acquire new values given by 

E^^E^ + {B^, - t^E.y , 4, = B^, + 2^^B,^phP . (184) 

but the gauge dependence of the corresponding contravectorial quantities 

E^^ = -f^'^E, , S'^ = \e>'''''B,p , (185) 

will be simpler, so much so that the vector B^ will actually be physically well defined 
in the sense of being independent of the choice of ether frame, since we shall have 

^ + ^I'^B^plf , B^^B". (186) 

The kinematic field 2-form closure condition 



V[;.F,,i = 0, (187) 

corresponds two of the 4-Maxwell equations, which are expressible in our 4-dimensionally 
covariant notation scheme as 



V/.S'* - , e^^'^'^W^Ep = -e^'V^B^' , (188) 

of which the second is interpretable as Faraday's law of magnetic induction. 

The trouble, in a Newtonian framework, is with the other two Maxwell equations, 
which specify the way an arbitrary source current 4- vector j^^ governs the dynamic 
evolution of the field. In a relativistic theory this is done by setting Vj^Fi^^ — 47rj^, 
where F^^^ is obtained from Fp_y by contraction with the non-degenerate contravariant 
spacetime metric g^^ . However the analogous Newtonian procedure of contraction with 
7^^^^ will, due to the degeneracy of the latter, give a result that is overdetermined, having 
a form that is expressible in terms of the rationalised magnetic field 

^ J_^MP^-ir^^ ^ J_^MP^-5^^ ^ (189) 
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which is ether gauge independent, and purely spacehke, 

H^"' = H^"' H^"'U = 0, (190) 

as the Ampere type equation 

V^H'"' = f . (191) 

It is evident from (190) that this can be satisfied only if the current is restricted to be 
similarly spacelike, in the sense of satisfying the consistency condition 

i% = 0. (192) 

The associated electromagnetic action density Ap, as similarly obtained from the usual 
relativistic action density F^^" F^^/lGn by substituting 'f^'^ for g^'^, will have the form 

Ap = -g'p , (193) 
where is the ether gauge independent magnetic energy density as given by 

g'p = , B'^ r),,B^B^ = \l''l'"'F,,F,, , (194) 

from which, by considering the effect of varying the degenerate metric 7^*^ for fixed F^^,^ 
the corresponding stress energy tensor is obtainable according to the ansatz (72) in the 
form 

T^,^H^PF,,-€J^. (195) 

The force density acting on the field (the opposite of the Faraday - Lorentz type elec- 
tromagnetic reaction on the relevant medium) will therefore be given by the expression 

/p. = V^T/, = fF^, . (196) 

This necessary restriction (192) is interpretable as meaning that there can be no net 
electric charge density, a condition which replaces the traditional Coulomb equation 
that would be expressed in our covariant notation scheme as V^E" = AnjH^, but 
which can be seen from (186) to be incompatible with Galilean invariance unless the 
magnetic part of the field is absent, and which even then will be incompatible with 
(184) except in the pure vacuum case for which there is no source current at all. 

To set up an ether frame invariant model for use as an exactly self consistent ap- 
proximation in a Newtonian framework we are thus faced with a choice between two 
generically incompatible alternatives. One possibility (which is likely to be most real- 
istic when insulating material is involved) is to use a scheme based on the Coulomb 
equation, which entails abandoning the Ampere equation and simply restricting the 
magnetic part of the field to be zero. The other possibility (more likely to be realistic 
for dealing with good conductors) which is the option chosen for the present work, is 
to use a scheme based on the Ampere equation (191) in conjunction with the force law 
(196), which entails abandoning the Coulomb equation and restricting the charge den- 
sity to be zero in accordance with (190). This effectively demotes the current from the 
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status of an independent dynamical variable to that of a derived quantity, and entails a 
concomitant loss of independence of the electric part of the field, which instead of the 
Coulomb equation, is required, in the most commonly used kind of model, to obey an 
Ohm type equation, of which the simplest version is covariantly expressible, in terms 
of the 4-velocity vector u^^ of the relevant supporting - fluid or solid - medium, as 
l^'^F^yU'^ = Kj'^, where k is a resistivity scalar that, in the case of a non-isotropic solid, 
might need to be replaced by a tensor. For positive resistivity, k, > 0, such an Ohm 
ansatz can be applied in the case of a composite medium involving entropy density as 
an independent degree of freedom, but its substitution in the force law (196) shows that 
it will entail a generically positive rate of energy transfer to the medium that will be 
given as a quadratic function of the (purely spacelike) current (191) by an expression 
of the form f^v = i^j'^, where = rj^yj^j'^. It is evident however that this will not 
in general be compatible with the identity (148) that must be satisfied for a single con- 
stituent medium of the kind to which the present study is restricted. To obtain a self 
consistent model involving just a simple solid or (barotropic) fluid, we need to restrict 
ourselves to the non-dissipative perfectly conducting case case for which the resistivity 
vanishes, ac = 0. 

10. Perfect magneto - elastic dynamics. 

It is evident from the foregoing considerations that the perfect conductivity con- 
dition needed to characterise a medium of the simple non dissipative kind considered 
here reduces to the perfect conductivity condition that is expressible covariantly as the 
condition 

F^.w'^ = 0, (197) 

which is interpretable as meaning that with respect to the local rest frame specified by 
the material 4-velocity the field is of a purely magnetic character. This condition 
is not just mathematically convenient but also justifiable - due to the relatively small 
mass of the electrons that are typically the main charge carriers - as a remarkably good 
approximation in many terrestrial applications and in a very wide range of astrophys- 
ical contexts, of which the most extensively studied so far have been those for which 
the the relevant material medium is a simple perfect fiuid, in which case the ensuing 
theory is what is known as perfect magnetohydrodynamics. As has been shown by the 
work of Jacob Bekenstein with Eleizer and Asaf Oron [27, 28], the elegant mathemat- 
ical properties of this kind of magnetohydrodynamic model arc easier to analyse in a 
fully rclativistic framework. Part of the motivation for the 4-dimensionally covariant 
approach developed here is to facilitate the extra work needed [27] for the treatment of 
the Newtonian limit. It is to be noted that the variational formulation developed below 
differs, in the fluid hmit, from the one developed by Bekenstein and Oron [27] who 
worked with Clebsch type potentials of the kind introduced in a relativistic context by 
Schutz [17]. The use of such Clebsch type variables is just one of several possibilities 
that may be convenient for various purposes in a purely fiuid context, but like most of 
the other alternatives it has the disadvantage of being unsuitable for generalisation to 
solids. For the purpose of setting up a variational formulation for the treatment of an 
elastic solid medium it has long been clear [15] that the only practical option is to work 
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in terms of world line displacements as characterised by comoving coordinate variable 
of the kind denoted here by . 

As a consequence of the closure property (187), it follows that the 2 - form will 
be "frozen in" in the sense of having vanishing Lie derivative with respect to the flow: 

u£ F^, = u'VpF^u + 2Fp^,V^]uP = . (198) 

since an antisymmetric matrix cannot have even rank, the orthogonality condition (197) 
implies that, as well as u^, the field F^^, possesses another independent null eigenvector, 
which can be taken to be as given by (185). A well known consequence is that the 
Maxwellian 2-form F^y will be conserved by Lie transport along any vector that is a 
linear combination of the form = aiB^ + a2U^ where ai and a2 are any scalar fields, 
and it can also be seen [16] that the 2-surface elements spanned by such vectors will 
mesh together to form a congruence of well defined fiux 2-surfaces. 

In the same way as remarked above about the stress tensor S^j^y, the world line or- 
thogonality property (197) is interpretable as meaning that F^^ naturally determines 
and is determined by a corresponding antisymmetric material base tensor with compo- 
nents F^g such that 

^/.. = ^AB?>f., (199) 

while the Lie transport condition (198) is interpretable as meaning that this induced 
field will be time independent. 

Fas = (200) 

SO that the covariant components will be those of a fixed 2-form field on the 3 
dimensional base manifold X . Morever it can be seen that as a consequence of the 
space-time closure property (187) the base space 2-form field will have a corresponding 
closure property, 

^Ub,c]=0, (201) 
which means that it will locally be expressible as the exterior derivative. 

Fab = 2A[,,,] , (202) 

of a "frozen in" 1-form field with components of the kind introduced in (90), whose 
spacetime pull back (88) thereby provides an expression of the familiar form 

F^, = 2V[^^,] . (203) 

It is to be remarked that this natural material gauge is not uniquely defined, since there 
is still some liberty in the choince of A^. As a result of the degeneracy property (197), 
there exists a current r]^ defined by 

T]" = l^'^'^^'M^F,,, (204) 

that is conserved in the sense of satisfying 

V^^M ^ = 0, (205) 
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and that has a time component 

V^t^, = -A^^B'' (206) 

which is proportional to the well known magnetic helicity scalar [29] . 

As a consequence of this "frozen in" behaviour, it can be seen that allowance for 
the effect of the magnetic field can be included directly within the perfect elasticity 
formalism developed in Section 5 simply by taking the energy density to have the form 

e^€o + €^ (207) 

in which € is the purely material part - depending just on the base coordinates and 
the induced metric components 7-^^ - that would remain when the field components 
F^g components are set to zero, while the other part is specified as a function not 
just of and 7-*^ but also of the field components F^^ whose status will be that of 
initial data that - subject to the closure condition (201) are freely specifiable, but that 
once chosen will evolve as fixed functions of the base coordinates q^. For a fixed base 
value of the base coordinates q^ the most general variation of & will determine not just a 
corresponding (symmetric) stress tensor 5'-*^ but also a corresponding (antisymmetric) 
magnetic field tensor by the prescription 

Se = - €^--)5^,, + \h--5F,^ , (208) 

which decomposes with 

gAB = gAB ^ gAB ^ (209) 

as the sum of parts given by 

56'o = |(5o""-€^o7.J<^7.., (210) 

and 

5€;, = i(5/- - €,Y'')51ab + kH-'^SF,, . (211) 

In the general case of a polarised medium the functional form of might be 
rather elaborate, but in the simple case of an unpolarised medium it is simply to be 
identified with the ordinary vacuum magnetic energy density as given by (194) which 
is translatable into terms material base fields as 

= \H--F,, , (212) 



where, consistently with (211) the magnetic field tensor H^'^ is given in this particular 
case simply by 

JJAB ^ l_^AC^BOp^^ (213) 

while the corresponding magnetic stress contribution will be obtainable from (211) as 

= H-cF^^ + €!,r, ■ (214) 
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which is equivalent to the expression given, using the notation (185), as 

= J-(25'^5''-5V)- (215) 

In order to consider the effect on wave propagation, as discussed in Section 8, we 
need the ensuing replacement of (177) for the tensor defined in the relativistic case by 
(172), which works out as 

= f 1 + ^^±p-] + -^{B^Y" - B^B-") . (216) 

We also need to evaluate the corresponding (unpolarised) magnetic elasticity contri- 
bution. This will be obtainable on the basis of the ansatz (162), which provides the 
expression 

J^ACBD — ^ACBD _j_ ^AB ^CD _ J^AB QCD _|_ JJDAJ^BC _|_ g< ^^AB^CD _|_ 2^'4[d^c]b^ 

_^25'jf [c^oJs _^ 25"/ [^7^1-* - _ (217) 

It immediately follows that the corresponding magnetic contribution in the character- 
istic matrix Q^'^ given by (175) will be expressible in terms of the relevant magnetic 
field tensor as 

Q^u ^ A^f'^'upu^ = -AniH^'PHp " + H^'PH'^'upU^) . (218) 

By introducing a spacelike unit vector that is chosen in the polarisation plane or- 
thogonal to (the unit vector in the direction of polarisation) in such a way that the 
magnetic induction vector will be expressible as 

B'' = ie^^^F.p = z/^ + E^t^ , (219) 

we can rewrite (218) as 

inQ^" = Bl^i'v'' - 2B^^B^u^^i''^ + B^^{j'"' - v^'v'') . (220) 

The complete characteristic matrix, 

QM- = Q^- + g (221) 

can easily be obtained in an explicit form if we suppose that the material contribution 
Qq^ has the form (180) that is relevant when the medium is in a simple isotropic state 
of the kind characterised by (178). In this (unpolarised, materially isotropic) case, we 
obtain an expression of the form 

Q^^- = Q^^ui^i,- - 2Qxi^(^i"^ + Q±(7''" - ^^^'') ■ (222) 

with 

«'=^° + T^ + §' ^x"^' «-=''° + i- P23) 
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The characteristic eigenvector equation (173) is easily soluble in the non-relativistic 
case, for which we simply have y^^^ = 7^*^: it can be seen that there will always be 
a transverse Alfven type mode, with polarisation covector ti^ orthogonal to both the 
propagation direction and the magnetic field direction, whose velocity v-,. will be given 

by 

^ p p Airp 

an expression that reduces to the well known formula v"^ = Hq/ p for propagation of 
transverse ( "wiggle" or "shake" ) modes in an isotropic solid when the magnetic field is 
absent. 

In the case of propagation parallel to the magnetic field direction, meaning oc 
so that Sj^ = 0, there will be a second (orthogonally polarised) transverse mode, with 
the same propagation speed v^, which in this case will be given by the expression 
■u^ = IJq/p + B^/Airp, which reduces to the well known Alfven formula = B^/Aup in 
the magnetohydrodynamic (purely fluid) limit in which the rigdity coefficient p vanishes. 
There will also be a purely longitudinal (sound type) mode with velocity v^^ that will 
be given, independently of the magnetic field strength, by the expression 

= ^ + (225) 

which reduces in the magnetohydrodynamic limit, — > 0, just to Newton's formula 
= dP/ dp for the speed of ordinary sound. 

More generally, when the propagation is not parallel to the magnetic field, the 
other two modes will be of mixed - partially longitudinal, partially transverse - type 
with polarisation in the plane generated by the propagation direction and the magnetic 
field direction, and with speeds and v_ that will be obtainable as the roots of the 
eigenvalue equation 

(g,i - pv^){Q^ - pv^) -Ql^O, (226) 

which gives 

2pvl = g,, + ± ^J{Q,^Q±y + ^Ql . (227) 

These solutions can be seen to be such that we shall have v+ — > v^^ and f _ — > v-^ in the 
limit of parallel propagation for which B± — * and B^^ — > S. 

In the relativistic case we shall still be able to use the same expressions (222) and 
(223) for Q'^" , but for the tensor y'^'^ in the characteristic equation (173) it will be 
necessary to use the less simple formula (216), which will be expressible, in a form 
analogous to (222), as 

pc^yf'-^ = - 2Q^i^^^'l''^ + Y^i^'i" + r(7''^ - uf'u'' - if'i") . (228) 

with 

l^,=cVo + §, Y^^c^p, + ^, Y^c^p, + ^, (229) 

where 

c% = cV + 6'o + Po. (230) 
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As before, there will always be a transverse Alfven type mode, with polarisation 
covector orthogonal to both the propagation direction and the the magnetic field 
direction, with velocity v-^ that will be given generically by 

vl_Q^_ 47r/.o + Bl 

In the case of propagation in the direction of the magnetic field, i.e. when B± — 
there will again be a second (orthogonally polarised) transverse mode, with the same 
propagation speed v^- , as well as a purely longitudinal (sound type) mode with velocity 
I'll that will be given in the generic case by the same formula as has long been well 
known [15] for the unmagnetised case, namely 

v: = l + ^. (232) 

For the generic case, in which the propagation is not parallel to the magnetic field, 
the relativistic generahsation of the equation (227) for the speeds v+ and v_ of the 
mixed modes - with polarisation in the plane generated by the propagation direction 
and the magnetic field direction - will be expressible in terms of the dimensionless ratios 

by 



l_Q, + Q±- 2QI ± ^( Q,|-Q^)2 + 4Q^(1-Q„)(1-Q7) 

2(1 ' ^ ^ 



r.2 on _ r)2 
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